Breakdown of Landau-Ginzburg- Wilson Scheme 
for Phase Transitions Between Orbital-Current Phases and Zero-gap Semiconductors 
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We show that Landau-Ginzburg- Wilson (LGW) scheme for symmetry-breaking transitions breaks 
down when spontaneous orbital currents cause quantum phase transitions from zero-gap semicon- 
ductors such as Dirac fermions to topological insulator (TI) or Chern insulator (CI). Topology change 
generates unconventional universality characterized by mean-field critical exponents j3 — n/ (d — n) 
and S = d/n arising from free-energy singularity |^| that induces enhanced quantum fluctua- 
tions. Here, d, and n are intensity of the orbital current, spatial dimension, and band-dispersion 
exponent, respectively. 

PACS numbers: 05.30.Rt,71.10.Fd,73.43.Lp,71.27.-l-a 



Introduction. — Critical phenomena of phase tran- 
sitions are classified into a small number of universal- 
ity classes. For such symmetry-breaking transitions, 
the LGW scheme is extremely successful Usual 
Fermi sea of electrons in crystalline solids, a vacuum 
for particle-hole excitations, is known to be unstable 
in the presence of electron-electron interactions, leading 
to spontaneous symmetry breakings. Competitions of 
different spontaneous symmetry breakings, for instance, 
magnetism and superconductivity are the subject of ex- 
tensive studies on strongly correlated electron systems. 
Although it has a rich diversity, most phase transitions 
accompanying the spontaneous symmetry breakings are 
primarily described by the LGW framework [l|, Q- In 
spite of its great success, however, the LGW scheme has 
recently been challenged from various viewpoints, espe- 
cially in strongly correlated electrons 

Contrary to the conventional Fermi sea, particle-hole 
excitations sometimes behave as massless Dirac fermions 
as realized in graphene and at surfaces of TI 
These zero-gap semiconductors attract much attention 
because of peculiar transport properties, such as ballistic 
transport 11 1 and quantum Hall effects [13, . Unlike 
the usual Fermi sea, band crossing points (Fermi points 
of the zero-gap semiconductors) are protected by certain 
symmetries j, 1411 and thus are inferred to have weaker 
instabilities ll|, ll5| . 

However, the zero-gap semiconductors are sensitive to 
the spin-orbit interaction, and time-reversal-symmetry 
broken fields such as magnetic flux. These stabilize topo- 
logically nontrivial phases as TI for the former and quan- 
tum Hall insulators (or CI) 



121 for the latter. It in- 



duces microscopic loop of spin (or charge) currents, for 
which we inclusively call orbital currents. We call phase 
transitions between the orbital-current phases (TI or CI) 
and zero-gap semiconductors orbital-current transition 
(OCT). 

Even without such external triggers, electron-electron 



interactions are proposed to cause spontaneous symme- 
try breakings generating the orbital currents 16^1^. 
This spontaneous-orbital-current phase is called a topo- 
logical Mott insulators (TMI). We call transitions to 
TMI, topological Mott transitions(TMT). 

In this Letter we propose that the universality class of 
OCT including TMT is unprecedented, because of the in- 
volvement of the topological change beyond the concept 
of the symmetry breaking, which causes non-LGW type 
singularities. This opens new avenues of studies on quan- 
tum critical phenomena and provides a clue for potential 
applications of TI (or CI) and its transitions. 

Differently from the Dirac fermions, quadratic band- 
crossing points have also been recently discussed [l7. 18. 
[20l . [2H . In general, zero-gap semiconductors may be de- 
scribed by a simple dispersion e_|_ and e_ as functions of 
momenta k, e± = ±u±fc", (n = 1 for Dirac fermions and 
n = 2 for quadratic band crossings) with a double degen- 
eracy at fc = 0. We present a comprehensive theory with 
dependences on n and spatial dimensionality d. 

In fact, Dirac fermions appear in honeycomb [l2| 
and diamond lattices [lij , while quadratic band-crossing 
points show up in kagome [13, liol , and pyrochlore lat- 
tices H, 2l|(FiglUa)-(c)), where, only in the pyrochlore 
lattice, three bands touch each other instead of two. 

The orbital currents make the band-crossing point 
massive (Fig[TJd)) in the simplest case as. 



(1) 



A mass m is scaled by, as we will discuss later, the orbital 
current multiplied by intersite Coulomb repulsion V; 
m (X VC. Orbital currents are sketched in Fig |ITe)- (g) 
(see also S.l of Supplemental Material for details (22|). 

Model and Method. — To understand the criticality of 
OCT, especially at T = 0, we expand the free energy F 
in terms of the order parameter C based on a microscopic 
hamiltonian and realistic band dispersions. Although we 
formally follow the spirit of LGW expansion, we will see 
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(a) honeycomb (b) Kagome (c) pyrochlore 




FIG. 1. (color online) (a) Honeycomb, (b) kagome and (c) 
pyrochlore lattice structures, (d) Schematics showing how 
degeneracy at a band-crossing point is lifted. The origin of 
the gap formation is understood from the lifting of the de- 
generacy related to the lattice symmetry. For example, in 2D 
lattices like the honeycomb and kagome lattices, the clock- 
and counterclock-wise rotating electrons are doubly degener- 
ate at the band-crossing point. The orbital currents break 
the symmetry, and lift the 2-fold degeneracy between them, 
leading to a gap. The loop current in this momentum space is 
equivalent to that in the real space shown in (g). Spin-orbital- 
current configurations on a hexagon (e) and a tetrahedron (f ) : 
Along each nearest neighbor bond {ri,rj), opposite spins on 
the directions ±6ij x dij flow in the opposite directions (g). 
Here, we define dij = ft — fj and bij = (fi + rj)/2 — Oij with 
Oij being the gravity center of the hexagon (e) (tetrahedron 
(f)). In (g), all the vectors schematically illustrate only the 
correct directions while the lengths are not correct. 



later that F has a smgularity expressed by that 
disrupts the LGW scheme. 

As a microscopic hamiltonian, we study an extended 
Hubbard model, 



(2) 



2>J 



where cj^. (Cj^.) is a creation (annihilation) operator for 
a cr-spin electron, = -I- n^^ is an electron density 
operator, is a pair of nearest-neighbor sites, and 

JJ (yjj) are an on-site (off-site) Coulomb repulsion. We 
only consider the nearest (Kj = Vi) and 2nd (Kj = V-i) 
neighbor interactions for Vij^ for simplicity. 

The Hartrce-Fock decoupling of the interaction (F) 
term in Eq.(I2]) is naturally described by a product of 



two spontaneous orbital currents that preserve lattice 
symmetries [l^ Among possible mean fields, two 

types of loop orbital currents exist: One is a mean field 
of the spin-orbital current. 



I fey X d^j\ 



(3) 



which creates time-rcversal-symmetric topological insu- 
lators and the other is a mean field of the charge-orbital 
current, 



2 X^jc'^jff 



bij X dij 
\bij X dij\ 



(4) 



which induces timc-rcvcrsal-symmctry broken Chcrn in- 



sulators. For definitions of bij and dij, sec Figlljg) (see 
also S.l of Supplemental Material 22]). We note that the 



form ^ is identical to the spin-orbit interaction uniquely 
allowed by the lattice symmetry and the time-reversal in- 



22|. 



variance 

Now in the Hartree-Fock theory, other types of orders 
such as the antiferromagnetic and charge-density-wave 
orders compete. In the plane of U, V and t, phase dia- 
grams clarified on various lattices after considering the 
competition indeed contain the region of stable symme- 
try breaking into cither Cs 7^ or (^c 7^ phases within 
realistic parameters [l6 . 
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Results. — The resultant free energy, /[(] is expanded 
by ( (for both C = Cc and ^ = C) as 

m - m = -K + + b±fs[C] + (higher order),(5) 

for small ^ close to OCT. Here a and b± arc constants. 
We add a spin-orbit coupling or magnetic flux A conju- 
gate to C, as a straightforward analogue of magnetic fields 
in magnetic phase transitions. In addition to the regular 
term proportional to a, the singular part fs[Q oc 
emerges with logarithmic corrections. For details of ex- 
pansion and derivation see Table |T] and S.2 of Supple- 
mental Material The coefficient for ^ > is not 
necessarily equal to 6_ for ( < 0. The expansion ^ is 
only piecewise analytic separately in ^ > and <^ < 
with a nonanalyticity at ^ = 0, originating from topo- 
logical nature of this transition, in contrast to analytic 
expansions assumed in the LGW framework. 

At nonzero temperatures T ^ 0, thermal fluctu- 
ations smear band-crossing points. The LGW form, 
f[C] - /[O] ~ -AC + AC^ + BC* + ©(C^) then recovers, 
leading to the Ising-type universality at T > 0. 

Now we describe generic features of the phase diagram 
for OCT obtained from Eq. ^ in Figll^a) and for the 
honeycomb lattice in Figl2jb). A semimetal (SM) is sta- 
bilized for a small interaction Va (a = 1,2) (correspond- 
ing to a > in Eq. ([5|)), at A = (along the white line in 
FigE^b)), while A drives continuous quantum phase tran- 
sitions to TI or CI at A 7^ 0. On the other hand, when 
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TABLE I. Singular part, order of transition, and critical exponents in doubly degenerate band-crossing point for several choices 
of band-dispersion exponent n and the spatial dimension d (for detailed derivation, see Supplemental Material S.2 First 
column shows fs, singular parts in f[(]. Second column shows order of transitions determined by following Ehrenfest's and 
Lifshitz's classifications. From 3rd to 5th columns, critical exponents around the QCP of TMT are shown. For (n, d) = (1,3), 

exponents (3 and 5 include logarithmic corrections denoted by ±0, where ~ ±^|a|/61n l/\a\ and ~ sign(A) \/|a[^ ^ In 1/|A| 
for \a\/t <ti 1 and |A|/t <^ 1. For (n, d) = (2, 2), "essential singular" means the instability by an infinitesimally small Va/t <^ 1. 
The pyrochlore belongs to (2,3), while it has some uniqueness partly because of the coexisting flat dispersions [2^ . We note 
that singularities share similarities to several antiferromagnetic transitions in zero-gap semiconductors [2^ . [2^ . 





fa[C] order of transition oc \a\ 
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FIG. 2. (color online) (a) Schematic phase diagram in the plane of a proportional to interaction Vc — V and conjugate field A. 
Free energy structures in each phase representing topological insulators (TI), (Chern insulators (CI)), TMI and semimetal (SM) 
are also shown, (b) Phase diagram of CI induced by Coulomb-repulsion on honeycomb lattices in the space of temperature 
T, A and interaction V2 obtained by Hartree-Fock calculation. For the honeycomb case, the Chern-type TMI (d 7^ 0, d = 0) 
and spin-orbit-type TMI (d 7^ 0, d = 0) are degenerate on the mean-field level, while the Chern-type will win because of the 
discrete symmetry breaking, against the SU{2) symmetry breaking of TMI, which is suppressed in 2D at T > 0. In the (pink) 
shaded area at A = 0, the time-reversal symmetry is spontaneously broken and its border at T 7^ (the solid (blue) line) 
illustrates its critical line. This critical line follows the Ising universality (or in other words, maps to the gas-liquid transition). 
The dashed (blue) line shows a crossover between the quantum and Ising regions characterized by <^ oc |a|, and ( oc lal'^/', 
respectively [23. The white circle at V — V2C — l-18t and A = T = is the unconventional quantum critical point (QCP) of 
the TMT from which the white topological critical line (semimetal) at T = and A = starts unlike the conventional QCP. 



Va exceeds the critical value Vac (for a < 0), sponta- 
neous orbital current (TMI phase) is stabilized for A = 
at T = 0, where the first-order transition separates A > 
and A < 0. The first-order jump terminates at a line of 
Landau's critical temperature Tc > This critical line 
(blue in Figl^t^b)) follows the Ising universality. 

Sandwiched by Landau's (blue) line at Tc > and 
the topological (white) line, a novel QCP (white cir- 
cle) emerges. The critical exponents of the QCP are 
generically /? = n/{d — n), 7 = 1 and 5 = d/n, in- 
dicating /3 > 1/2 and (5 < 3 (see Table |T| , completely 
opposite to the standard predictions of LGW symmetry- 



breaking transitions, which always satisfy /3 < 1/2 and 
(5 > 3. Here, /3 and 5 are basic critical exponents 
to characterize how the order grows, as Q oc \a\^ at 
A = 0,a < and ( cx |A|^/'' at a = 0, respectively. The 
present unusual universality is governed by the topol- 
ogy, similarly to the purely topological ones identified as 
the marginal quantum criticalities for Mott and Lifshitz 
transitionsjllilliij. 

Strong correlation effects emerge as suppressions of 
simple orders with residual entropy and large quantum 
fluctuations at low temperatures. In the conventional 
systems with large Fermi surfaces, they emerge typically 
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as momentum and orbital differentiations as found in 
pseudogap and Fermi arc formation in the cuprate su- 
perconductors [26, 27 1 . In the zero-gap semiconductors, 
instead, large quantum fluctuations emerge quite differ- 
ently as the unusual exponents such as large /? and small 
S, resulting in a slow and suppressed growth of the or- 
der, although the transition is strictly protected by the 
topological distinction at T = 0. 

Now we consider an exceptional but interesting exam- 
ple of the pyrochlore lattice. Although the pyrochlore 
lattice has a unique feature [13], /[C] follows the same 
form as Eq. © with = \Cs\^/'^ and < 6+ < 6_. 

The obtained phase diagram in Figl3]contains several dis- 
tinct phases; SM, TMI (at A = 0, > Vic), TI (at A > 0) 
and TMI coexiting with charge order (TMI-hCDW) (at 
A < 0,T = 0). AU the phases meet at the QCP (white 
circle), which follows a novel universality similar to the 
other lattices (see Tabic |l| . The extension of a critical 
line at T = from the QCP (white lines in FigsIS] and 
[3]) is common protected by the topological nature [5[ and 
is fundamentally different from the LGW-type criticality 
even in itinerant electron systems j28l - l30j that does not 
have such a critical line and is simply described hy d + z 
dimensional classical LGW transitions, where z is the dy- 
namical exponent. An important difference of Figl3]from 
the symmetric case in Fig[5] is that the SM persists for 
A < 0. In addition, a first-order transition surface (light 
green surface) (represented by T = Tist(A, Vi)) bends to 
A > 0, which separates the "topological semiconductor" 
phase at r > and A > into two: Liquid-like {T < Tist) 
and gas-like (T > Tist) semiconductors. 

It is intriguing to associate the first-order transitions 
between gas- and liquid-like topological semiconductors 
(for A > 0) as well as between SM and TMI+CDW (for 
A < 0), to puzzling metal-semiconductor transitions in 
many pyrochlore compounds with bad metallic behavior 
as found in TI2RU2O7 [I8, 31| and Ln2lr207[3l| for sev- 
eral lanthanoid elements Ln. Here, we have ignored the 
orbital degeneracy and orbital dependent anisotropy of 
transfers of the d bands, while in many cases basic struc- 
tures at the band crossing are preserved even with this 
complexity. Indeed, triply degenerate bands at F point 
near the Fermi level shown in the local density approxi- 
mation [sij supports the relevance of the present general 
theory. Whenever zero-gap semiconductors are found, 
our general scheme for (n, d) applies. 

A nonzero spin-orbital current |C| 7^ in the TMI 
breaks the rotational SU(2) invariance of the global spin 
quantization axis preserved in the original Hamiltonian. 
Therefore, the Nambu-Goldstone mode of this order is 
the spin rotation coupled to the orbital motion. We call 
this new collective excitation, spin-orbiton. The exis- 
tence of the spin-orbiton is a way to distinguish TMIs 
from simple TIs induced by the spin-orbit couplings. De- 
tecting whether such an excitation exists is an experimen- 
tal challenge. In the SM phase, diverging spin-orbiton 
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FIG. 3. (color) Phase diagram of pyrochlore lattice. When 
A = 0, the system undergoes a transition from SM for small 
Vi to TMI (dark orange line) at the QCP, Vi = Vic ~ 2.62t 
at T = (white circle). The QCP shows a new universality 
similar to Fig[2] beyond the LGW scheme. When A < 0, 
SMs for small Vi make a first-order transition (the light 
blue surface), to a topological semiconductor coexisting with 
charge order (TS-fCDW). The TMI coexiting with charge or- 
der (TMI+CDW) at A < 0, r = crossovers to TS+CDW at 
finite temperatures. It terminates at a quantum critical line 
(dark blue). For A > 0, a topological insulator with a gap is 
stabilized for all Vi > at T = 0. At T > 0, through a (light 
green) surface, ^ jumps. This surface separates a gas-like and 
liquid-like topological semiconductors (TS) and terminates on 
the Ising critical (bold green) line, which further terminates 
at the QCP. The QCP also terminates two quantum critical 
lines at T = 0, one along \ = 0,V < Vic (bold white line) 
and the other (the dark blue line) representing that between 
TS-f CDW and SM. Here we note that, in the 3D lattices as 
the pyrochlore, the Chern-type TMI is not allowed and only 
the spin-current phase with ("s 7^ is allowed [l^] . 



fluctuations may couple to particle-hole excitations re- 
sulting in an overdamped modCjand alter the critical 
dynamics (dynamical exponent) [29|. This is an issue to 
be pursued. 

Furthermore, spin-orbital currents that preserve the 
time reversal and lattice symmetries, are invisible in 
many standard experiments. TMI thus is a potential 
candidate of "hidden" orders behind mysterious insulat- 
ing behaviors, as the orbital-current-induced pseudogap 
scenario discussed in cuprate physics [33|, |34 1 . 

Summary. — Our findings present a solid foundation 
for unconventional critical phenomena of OCT, and serve 
a basis for utilization of the OCT such as switching be- 
tween the topologically distinct phases. For more quan- 
titative estimate of the critical exponents by consider- 
ing fluctuations, for example, fermionic renormalization 
group [3^ or a renormalization group method applicable 
to non-analytic free-energy expansions will be helpful and 
are left for future studies. Strong quantum fluctuations 
expected around the present OCT may offer a basis for 
unprecedented quantum phases mediated by the spin- 
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orbitons. 
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Supplemental Material 

for "Breakdown of Landau-Ginzburg- Wilson 
scheme for phase transitions between orbital- 
current phases and zero-gap semiconductors " 
by Moyuru Kurita, Youhci Yamaji, and Masatoshi Imada 

In this Supplemental Material, we show theoretical 
methods used to derive our results in the main article. 
First, we show that the form of the orbital-current mean 
field is uniquely determined if we preserve the lattice 
and time-reversal symmetries. We next derive the free 
energy expansions for general degenerate band crossings 
and apply them to some specific examples, a honeycomb, 
kagome and pyrochlore lattices. As is introduced in the 
main article, we start from an extended Hubbard model 
on these lattices, 

(S.l) 



S.l ORBITAL-CURRENT MEAN FIELD 



(a) Honeycomb (b) Kagome (c) Pyrochlore 




(d) , , (e) (f) 




Fig. S 1. (a)-(c) Examples for vectors fi, fj, Oij, bij, and 
dij for the honeycomb, kagome and pyrochlore lattice, (d)-(f) 
Examples of symmetric operations applied to a bond 



Here, we show that the spin-orbit interaction as well 
as the orbital-current mean-field has a unique form if we 
preserve the lattice and time reversal symmetries. Let 
us first consider some symmetry operations satisfied by 
honeycomb, kagome, and pyrochlore lattices. When we 
pick up a bond between i-th (f^) and j-th sites (fj) of a 
lattice, the system is invariant under 3 important sym- 
metric operations related to the bond: (1) The time re- 
versal operation O, (2) a refiection with respect to a plane 
perpendicular to the bond and containing the center of 
the bond, which is denoted by i?^ (3) a 180°-rotation 
around the axis perpendicular to the bond and crossing 
at the center of the bond each other, which is denoted as 
°s . For all the lattices considered, the axis for i?}f"^ 

that keeps the lattice invariant is along the vector bij 
pointing to the center of the bond from the gravity 
center Oij of the unit cell containing i and j (see Fig.^. 
This is defined as bij = {n + rj)/2 — oij. For the honey- 
comb, kagome and pyrochlore lattices, Oij is the gravity 
center of the hexagon, triangle and tetrahedron, respec- 
tively, as illustrated in Fig.^a)-(c). 

Next we introduce a current fiowing through the bond 



along 



rj, described by the hermitian 

f(a) It 



Let us examine symmetric properties of "/(f])- As 
charge currents break the time-reversal symmetry while 
spin currents preserve it, J,^°\ and J^i"'.^'^' are time- 
reversal odd and even, respectively. 

Similarly, the reflection R'^ and the rotation -Rjf^-) 

transform into -,[M,] = Rm°[Ms] = 

— J(^^jy Therefore we note that, if the lattice symmetry 
satisfies these refiection and rotation as symmetric oper- 
ations, the charge-current mean field breaks the reflec- 
tion/rotation symmetry in addition to the time- reversal 
symmetry. When the charge currents form a loop, after 
setting currents on every bond as we will discuss below, 
these currents break the so-called chiral symmetry. 

On the other hand, spin current operators are 
transformed in complicated ways: The reflection i?^ --^ 
makes a spin rotate 180° around dy-, while -R(f°) makes 

the spin rotate 180° around bij. In addition, the both 
operations exchange the sites i and j. For example, if 



(-|dy |,0,0) and bij = (0, |6y |,0), the 
"i^) fiv) ji^) ^ 



we take dij - 

reflection R'^i j-j transforms a vector (J^ 



charge/spin current operator ["^(i^j)]^) defined by into {-jlij)^J(ij):Jlij))- while the rotation i?(,f°) trans- 



4"' 



• (S.2) 



Here, a 2x 2 matrix ct" represents, for a = 0, a 2D identity 
matrix describing the charge current or, for a = x,y^z^ 
Pauli matrices describing the spin currents. 



forms it into {jjj^^jy ^^^'^^ ^ sp^n current 

preserves the two symmetries, the spin alignment vector 
of the current from j-th to i-th sites, has to be 
along the z-axis. More generally, when the above three 
symmetries are preserved, the unit vector j) has to be 
perpendicular to both of dij and bij and uniquely has the 
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form n 



±btj X dij/\bij X di. 



Now, the spin current operator is expressed by 



(S.3) 



We also note that for the honeycomb lattice, the bond 
has to be between the 2nd neighbor sites, because 
the nearest-neighbor and the 3rd-neighbor bonds have an 
additional 180°-rotation symmetry with the axis crossing 
the honeycomb plane perpendicularly at the center of the 
bond. This symmetry prohibits the currents even in the 
bij X dij direction. In a similar way, for the kagome and 
pyrochlorc lattices, dij ^ (i, j) has to be between the 
nearest-neighbor sites. (See Figllja)-(c).) 

After choosing the direction and spin alignment of the 
flow, we apply a symmetry operation, which belongs to 
the point group of the underlying lattice, to the bond 
For example, in the honeycomb lattice, every 2nd- 
neighbor bond, where orbital currents flow, is involved 
in a hexagon. A 6- fold rotation (Cg rotation), around 
the center of such a hexagon, Oij, is one of the examples 
of symmetric operations. If we apply the Cg rotation, 
we obtain the bond shown in Fig[TJd). After suc- 

cessive applications of Cg rotations around Oy, we can 
determine the direction and spin alignment of the flow of 
5 more bonds involved in the hexagon. Under these op- 
erations, dij is invariant, and bij and dij are transformed 
into, for example, bi'ji and diij' , respectively. If we fur- 
ther use the translational invariance of the honeycomb 
lattice, we can determine the direction and spin align- 
ment of the flow for all of the 2nd neighbor bonds on 
the lattice. For example, we obtain the bond by 
translating the bond j) shown in Figlljd). Simultane- 
ously, vectors bij and dij are transformed into bi"j" and 
di-j" as well. When the unit cells share their corner, for 
a kagome and pyrochlorc lattice, the inversion operations 
shown in FigHJe) and (f) can be used. 

Then an explicit expression for orbital-current mean 
fields with a magnitude ^, together with a hopping rcnor- 
malization g due to a usual Fock term is obtained from 
Eq. (jS.2l) for the spin-orbital-currcnt mean field 



\bij X dy l 



(S.4) 



which creates the topological insulators, or a charge- 
orbital-currcnt mean field. 



= .9 - < 



bij X dij 

\bij X d,j\ 



Sap, (S.5) 



which induces the Chern insulators (absent in 3D lat- 
tices). These mean fields arc generated through a decou- 



pling of the two-body interaction term. 



(S.6) 



Here we note that the orbital-currents instabilities 
compete with charge and/or spin density wave instabil- 
ities in real lattice structures. Many studies were car- 
ried out to search the conditions that an orbital-current 
phase becomes a stable mean-field ground state as dis- 
cussed by S. Raghu, X.-L. Qi, C. Honerkamp, and S.-C. 
Zhang, Phys. Rev. Lett. 100, 156401 (2008), J. Wen, 
A. Riiegg, C.-C. J. Wang and G. A. Phys. Rev. B 82, 
075125 (2010), and M. Kurita, Y. Yamaji, and M. Imada, 
J. Phys. Soc. Jpn. 80, 044708 (2011) (Refs.16-18 of main 
article). For instance, the spin-orbital currents are sta- 
bilized for reasonably realistic values, 4.5 ^ U/t <6 and 
2.5 < Vi/t < 3 for a pyrochlorc lattice with one electron 
per site. 



S.2 FREE ENERGY EXPANSION 

Based on orbital-current mean fields given in Eqs. 
(|S.4p and (|S.5p . the mean-field hamiltonian for specific 
lattices is derived, with an appropriate choice of param- 
eters, for example, ([/, t^i, V2). First, we derive a general 
and simple analytic formulae of free-energy expansions 
for doubly-degenerate band crossings, by focusing on low- 
energy band dispersion around band-crossing points. It 
is applied to the case of the honeycomb and kagome lat- 
tices. Second, we give a free energy expansion for a triply- 
degenerate band crossing and applies to the pyrochlorc 
lattice. We also refer to a peculiar inversion of valence 
and conduction bands in a pyrochlorc lattice occurring 
by the sign change of the spin orbital current, Here, 
we note that the transition discussed in this article is 
different from that between the topological and band in- 
sulators discussed by S. Murakami, New J. Phys. 12, 
065008 (2007). 



Doubly degenerate band crossings 

When magnitude of an orbital current C. between a-th 
neighbor sites becomes finite, two bands split in the mo- 
mentum space, e±[C]{p), with a gap scaled by Then 
we obtain a mean-field free energy per unit cell as. 



m 



2kBT 



El 



In 



1 + e" 



^(e+[C](p)-M)/fcBT 



In 



1 + e 



-ie.lC]{p)-tJ.)/kBT 



(S.7) 
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where Nu is the number of unit cells, Za is the coordina- 
tion number of the a-th neighbor site, is the number 
of sites in an unit cell, and Va is the Coulomb repulsion 
between a-th neighbor sites. The chemical potential ^ is 
chosen to keep electron density or filling listed in Table ^ 
to keep the valence band e_ fully-filled at T = 0. 

Here we focus on a small energy window around a band 
crossing point pbc, and introduce a cutoff A, which does 
not change any essential physics. If we choose a proper 
A for each lattice, full numerical results for the original 
microscopic hamiltonian are quantitatively reproduced. 
Then we can expand the band dispersion e± with re- 
spect to the momentum measured from the band crossing 
point, k as 



e±[C]{pBC + k) ^ eik) ± 



vl\k\ 



2n 



(S.8) 



where a "mass" m is given as m — c^VaC with a coeffi- 
cient Cq. The dispersion e(A:), coefficient Cc, and "veloci- 
ties" v± depending on the lattices are listed in Table ^ 
At r = 0, the free energy has a simple analytic form. 



F[C] - F[0] 2nBc^i 



Nu 



Vb2 



Zcn^Vc^C. (S.9) 



Q 



k'^-^dk 



(S.IO) 



where k = |fc|, nee is the number of band-crossing points, 
Fbz is the volume of the Brillouin zone, and = Stt (47r) 
for d = 2 (d = 3). 

Singular part of the free energy arises from 







1 


1 - 


i" 


2 9t? 



/•A/|mh 



(S.ll) 

where we rescale the integration variable as g = /c/jmp/" 
to absorb the singular m dependence out of the inte- 
grand. Here we note that the integral in the right hand 
side of Eq. (jS.ll[) causes ultra-violet divergences in the 
limit |to| — > 0-I-, for n < d, which is compensated by the 
prefactor |to|'^/"^^ but leaves the singularity. In fact, by 



expanding the integrand in power of l/g^", Q turns out 
to have series of singular terms oc |to|2^+^~'^/" with any 
integer ^ > 0, which yields a singular expansion in terms 
ofC. 

For example, for the honeycomb lattice with (n, d) — 
(1,2), namely, for d/n+1 ~ 3, we obtain a non-Ginzburg- 
Landau- Wilson-type free energy expansion as 



F[C] - F[0] 



aC + b\C\^ + (higher order), (S.12) 



where the coefficients a and b are given in Table ^ 

For d = {2i+ l)n, the integral Eq. ((Sll|) yields a 
logarithmic correction proportional to ln|C|. Indeed, for 
n ~ d, we obtain the free energy expansion as 



F[C] - F[0] 
Nu 



6C^ In Id + aC^ + (higher order). 



(S.13) 



where the coefficients a and b are given in Table ^ for 
the kagome lattice with (n, d) = (2, 2). 

The band crossing satisfying d = 3ri or (n, d) = (1, 3) 
is not known in the real lattice to our knowledge. We 
can, however, predict a general form of the expansion as 



F[C] - F[0] _ 2nBz^d 
Nu 



Vbz 4du 



Vbz 



(higher order). 



c2 yV2 

Q *^ a S 



(S.14) 



Triply degenerate band crossings: pyrochlores 

Now we give a free-energy expansion for the pyrochlore 
lattice, as an example of the case (n, d) = (2, 3). 

The mean-field hamiltonian for topological Mott states 
is given as 



4-^ 



where the 8x8 Bloch hamiltonian h o{k) for a non- 
interacting electrons on a pyrochlore lattice is given by 



Vo(fc) 



cos(fca; — ky) cos(fcj, — kz) cos{kz + kx) 

cos{kx — ky) cos{kz — kx) cos(fcy -I- k^) 

cos(fcy — kz) cos{kz — kx) cos(fc^ -I- ky) 

cos{kz + kx) cos{ky + kz) cos{kx + ky) 



(S.16) 



with a spin-independent renormalized hopping matrix t' = t + Vi (cj^Cj^'^ and a 2 x 2 unit matrix I2 acting on the 
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Table S I. Summary of band dispersions, parameters depending on lattice structures, and coefficients for the free energy 
expansions. The cutoff A for the kagome lattice {{n,d) = (2,2)) is numerically not accurately determined yet, because of a 
numerical difficulty in the essential singular behavior. 



filliuE 



e(fe) 



a 
Vbz 



a 

b 
A 



honeycomb (n, d) = (f,2) 



kagome {n, d) — (2, 2) 



one electron per site 



4/3 electrons per site 




3^/3 



2 
6 
2 

2 



12V2 



b± 



0.590 



2t' - 
2^/3 

t'/2 



12Vi 



b± 



1 
4 
1 
3 



l + 21n(^ 



pyrochlore (n, d) — (2, 3) 



one electron per site 



241/1 - 2^ 



9.37t' 



2 6_ = 55At' 
0.507 



spin degrees of freedom. The spin orbital current j ^(fc) is given by a 8 x 8-matrix, 






'^xy 


Uyz 


Vzx 







Uzx 


Vyz 




-Uzx 





'^xy 


-Vzx 


-Vyz 


Vxy 






with 



""tj 

Vii 



(o-j +(Tj)c0s(fei 
: ((Ti - (Tj) COs(fc, 



(S.17) 



(S.18) 
(S.19) 



Here $t 



(c - ,c - ,c - ,c - ,c - ,c - ,c - ,c - ) represents a veetor notation of the ereation operators of the 

^ Ifet Ifci' 2fct 2fc4.' 3fct SfeJ.' 4fct 4A:i' ^ 

Bloeh states, with the suffix 1 to 4 in the creation operators representing the 4 sites of the unit cell. 

I 



By diagonalizing the mean-field hamiltonian, we ob- 
tain 4 bands. Out of the 4 bands, we focus on 3 bands 
touching each other at the F point. In contrast to sim- 
ple conduction- and valence-band dispersions already dis- 
cussed, here, we need to handle 3 bands for C = 0, 



e = 2t',2t',2t'-2t'fc2 + e)(fc4 



(S.20) 



By considering the spin degenracy, we have 6 bands 
touching at the band crossing point. 

In the presence of non-zero spin orbital currents C 7^ 0; 
we obtain 3 split bands. 



e = 2t' + 2^/21/lC, 2^' - V2ViC - t'k^ 



± 



y^(V2FiC-t'fc2)2 



(S.21) 



where the hopping matrix t is renormalized as t' = t + 
Vi (cj^Cj^'^. Each dispersion is doubly degenerate. At 
the r point of the Brillouin zone, fc = 0, in terms of 
the group theory, 6 degenerate bands in the T ® SU{2)- 



manifold split into a doublet (-E5/2) and quartet 
due to spin orbital currents. 

As it has been already mentioned by H.-M. Guo, and 
M. Franz, Phys. Rev. Lett. 103, 206805 (2009) (Ref.21 
of main article), for non- interacting electrons with spin- 
orbit couplings on pyrochlores, the sign of C is crucial 
whether the system becomes topological insulators or 
not. We note that, for the honeycomb and kagome lat- 
tices, topological Mott states appear for both positive 
and negative C. However, for the pyrochlore lattice, the 
sign of C. reverses the level of the doublet /2 a-nd quar- 
tet G3/2- Only if the doublet i?5/2 becomes occupied, the 
system becomes a topological Mott insulator. 

Then the free energy at T = is given by 

.A , 



m ~ m 



32 



24^1 C^. 



(S.22) 
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From this form of the free energy, we have 



■(C = 0) = -^1^ + 48^1 (S.23) 



and 



1536i' fV^ 



x'^(18- sgn(C)y2a;2) 



dx 



(18-sgn(C)2\/2x2 +a;4)i 



(S.24) 



from which the energy expansion is given by 

^ - ^ + ae + b^\C\i + 0{C% (S.25) 
where a and b± are summarized in Table ^ 



